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5.1  THE  IMPULSE  APPROXIMATION 

Following  the  historically  significant  papers  of  DuMond.1,2  almost  all  Compton 
profile  experiments  have  been  analysed  in  the  so-called  ‘impulse  approximation’.  It 
is  only  in  this  approximation  that  a simple  relationship  exists  between  the  Compton 
profile  J(q)  and  the  electron  momentum  density  (EMD)  of  the  bound  state  of  the 
atom,  molecule,  or  solid  being  probed.  Recent  review  articles3,4  have  derived  the 
impulse  approximation  from  the  Waller  Hartree  theory3  or  from  the  time  dependent 
operator  formalism  results  first  obtained  by  Eisenberger  and  Plat/man6  (see  chapter  2). 
The  authors  of  previous  review  articles  seem  to  feel  that  ‘while  corrections  to  the 
impulse  approximation  can  now  be  made,  they  tend  to  be  small'4  and  that  a 
‘laborious  calculation  of  the  final  state  wave  function  is  unlikely  to  be  necessary 
except  for  deeply  bound  core  electrons,  or  for  very  accurate  experimental  data’.3 
It  is  the  purpose  of  this  first  section  to  present  a more  straightforward  derivation 
of  the  impulse  approximation  and  to  comment  in  more  detail  on  the  accuracy  of  this 
approximation. 


5.1.1  Derivation  of  the  impulse  approximation 

In  the  non-relativistic  Schrodinger  equation,  the  perturbation  on  an  electron  due  to 
the  presence  of  a magnetic  field  w ith  vector  potential  A can  be  w ritten 

e e2 

H' = p.A+- j A.  A (5.1) 

m0c  2m0c~ 

where  p is  the  electron  momentum,  e is  the  electron  charge,  m0  is  the  electron  rest 
mass,  and  c is  the  speed  of  light.  In  a scattering  process,  the  photon  must  suffer 
annihilation  and  then  creation.  By  expanding  A in  creation  and  annihilation 
operators,  it  is  clear  that  the  scattering  process  is  of  order  A2.  Thus,  for  a first  order 
perturbation  calculation  the  p . A term  contributes  nothing.  It  is  only  in  second  order 
that  this  term  can  contribute  to  order  A2.  However,  if  one  assumes  that  the  incident 
photon  energy  lies  significantly  above  the  ionization  edge,  the  energy  denominators 
in  the  second  order  perturbation  theory  expressions  will  make  such  p . A contributions 
quite  small.  Some  estimates  of  p . A contributions  have  been  made  for  the  hydrogenic 
problem.'-8  One  usually  considers  only  the  A2  term  in  first-order  perturbation 
theory  in  the  calculation  of  Compton  scattering.  The  scattering  process  can  be 
compared  with  the  absorptive  photo-electric  effect  where  only  the  destruction  of  a 
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photon  is  required  and  only  the  p . A term  is  considered  in  lowest  order  perturbation 
theory  (first  order  in  A). 

Using  Fermi's  golden  rule  for  calculating  the  transition  probability  and  cross- 
section  along  the  lines  of  Schwinger*1  leads  to  the  cross-section  cited  by  Eisenberger 
and  Plat/man6  for  a single  eleetron  in  an  independent  electron  model 


where  k,  = (Oi/hc  and  <p  is  the  scattering  angle.  ipt  and  i//2  are  the  wave  lunction; 
for  the  electron  in  its  initial  bound  state  and  its  final  continuum  state  respectively 
The  Thomson  cross-section  for  unpolari/ed  incident  radiation  is  given  by 


Eisenberger  and  Plat/man6  include  a factor  (aii/wi)2  in  their  definition  of  (drr , dQ)r 
which  modifies  their  eq.  6 to  include  the  factor  aq/roj  rather  than  wj/wi  in  our 
eq.  (5.2)  above.  The  delta  function  in  eq.  (5.2)  guarantees  conservation  of  energy.  In 
the  general  scattering  case  no  conservation  of  momentum  between  the  photon  and 
electron  is  required  as  the  nucleus  itself  may  carry  away  momentum  during  the 
collision. 

To  derive  the  impulse  result  let  us  assume  in  the  f)  function  that  the  electron 
can  be  treated  as  free  (but  moving)  during  the  entire  collision  process.  Before  the 
collision  let  the  electron  have  momentum  p,  and  after  the  collision  momentum  p>. 
Then 


Conservation  of  momentum  for  this  free  electron  case  gives 

p2  = Pl  + /ik;  pi  ~ p\  + 2/ik  . p,  + h2k 
and  the  d function  becomes  (see  eq.  (2.4) ) 


Within  the  impulse  approximation  ip 2 is  taken  to  be  a plane  wave 

i/f2(r)  =■  eip'-,!  *.  (5.8) 

The  momentum  space  bound  state  wave  function  / 1 i p i ) is  the  Fourier  transform  of 
the  corresponding  position  space  function  so  that 
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X,(Pi)  = (-n/i)  •' : Ji//,(rlc  ■' * dV  (5.9| 

and  the  sum  over  final  states  is 

£-*  |(2afc)-3d3/»2.  (510) 

Choosing  k to  lie  in  the  r direction  and  using  momentum  conservation 


X|<i/'i|e'k-r|V'2>!2^(£2  - Ei  - a>) 


= fd3p2|(2,tA)-«|d3re*- '^(r) |2  J [p,,  - 


Defining 


eq.  (5.1 1)  reduces  to 


<?  = 


imv  hk 

hk  ~ ~2 


dV2|Zi(Pi)|2<5^  (Pis  - q)  )• 


(5.1  II 


(5.12) 


(5.13) 


Noting  that  for  to,,  (p,  and  co2  (or  tu)  fixed,  k is  fixed  as  well,  it  follows  from  eq.  (5.6) 
that  the  integral  over  d 3p2  is  equivalent  to  an  integration  over  d3p,  dplA  dp,,  d 
Making  use  of  the  properties  of  the  <5  function,  eq.  (5.13)  becomes 

+ OO 

\Xi(Pix.Piy.q)  |2  dp u dp,,.  (5.14) 


which  is  equivalent  to  an  integration  over  the  plane  in  momentum  space  p,;  = q. 
For  a momentum  density  which  is  spherically  symmetric,  it  is  convenient  to  rewrite 
eq.  (5.14)  in  cylindrical  coordinates  from  which  it  reduces  to  (miiik)J(q)  where 

J{q)  = 2ti  j |/i(p,)|2p,  dp,.  (5.15) 

Jf«l 


Fxpressing  p and  q in  atomic  units  of  momentum  h/o0.  eq.  (5.2)  becomes  the  well 
known  result  for  the  Compton  profile  in  the  impulse  approximation 


d:<7  /dff\  W;  IIUIq  1 

dQ  dm  \dfi/Tcu,  /i2  ka0  ^ 


(516) 


The  derivation  has  been  presented  in  this  much  detail  to  bring  out  clearly  all  the 
approximations  leading  to  the  impulse  result  of  eqs.  (5.15)  and  (5.16).  and  also  to 
serve  as  an  introduction  to  the  units  commonlv  used  in  Compton  work. 

From  the  foregoing  analysis  it  should  bequite  clear  that  the  impulse  approxima- 
tion is  equivalent  to  a photon  scattering  inelastically  from  a free  electron  gas  with 


Ain  u.  I0> 


a spherically  symmetric  (although  this  is  not  intrinsic  to  impulse)  momentum 
distribution  where  both  energy  and  momentum  are  conserved  during  the  collision. 
The  momentum  density  for  any  particular  momentum  in  the  free  electron  gas  is 
obtained  from  the  square  of  the  Fourier  transform  of  the  initial  bound  state  radial 
wave  function.  The  energy  momentum  conservation  relations  determine  that  for  an 
incident  photon  with  energy  o> i scattering  at  any  fixed  k (r/>  and  ei2),  contributions 
to  the  scattering  from  electrons  with  momenta  pt  can  occur  only  if  the  projection  of 
pt  on  the  scattering  vector  k equals  q.  This  can  only  occur  if  /q  is  greater  than  or 
equal  to  q. 

DuMond1  intuitively  represented  the  atom  as  a free  electron  gas.  Solving  the 
energy-momentum  conservation  equations  relativistically,  he  was  able  to  show  that 
because  the  electron  moves  with  momentum  p,  before  the  scattering,  the  wave- 
length of  the  scattered  photon  is  Doppler  shifted  from  the  usual  Compton  result 
according  to 


A 2 — At  = - (1  cos  <t>)  + q/.*jmnc 

m0c 


(5.17) 


Where 

A*  (A i 4 A2  - 2AiA2  cos  (/>) 1 2 (5.18) 

and  q is  the  projection  of  the  initial  momentum  on  the  scattering  vector  as  before. 
This  result  occurs  after  setting  the  initial  electron  energy  to  m0c2  and  may  not  he 
immediately  discerned  from  DuMond’s  papers  since  he  defines  A*  and  q slightly 
differently  from  the  above.  Since  the  A*  term  varies  slowly  over  the  profile,  eq.  (5.17) 
shows  that  A2  varies  linearly  with  q.  Equation  (5.17)  suggests  that  a large  and 
negative  q corresponds  to  a small  A 2 or  a large  m2  of  the  scattered  photon,  whereas 
a large  and  positive  q corresponds  to  a large  As  or  a small  an  of  the  outgoing 
photon,  and  a corresponding  high  energy  for  the  ejected  electron.  From  eq.  (5.17), 
one  can  define  q as 

q = m0c!//.*  (5.19) 

where 

/ = A2  - A]  — (/i,'m0c)(1  - cos  (j> ) (5.20) 

represents  the  displacement  in  wavelength  from  Compton  scattering  from  a free 
electron  at  rest  In  his  anal>sis  DuMond  first  considered  the  scattering  from  a 
mono-energetic  isotropic  free  electron  gas  with  momentum  p{.  It  is  clear  from 
eq.  (5  17)  that  such  a gas  would  give  rise  to  a flat  intensity  distribution  with  a width 
in  )2  of  Ipx/S/mc.  It  also  follows  that  the  area  of  this  intensity  band  is  proportional 
to  the  number  of  electrons  with  momentum  pt.  Then  considering  a gas  with  a 
distribution  in  momentum  and  assuming  scattering  equally  probable  from  any 
momentum,  DuMond  obtained  the  result  of  eq.  (5. 1 5)  for  the  unnormali/ed  intensity 
distribution.  Of  course  this  omits  the  slowly  varying  factors  c>:'nilk  contained  in  the 
correct  impulse  result  given  in  eq.  (5.16).  It  is  certainly  a tribute  to  DuMond  that  his 
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intuitive  approach  brought  out  the  most  salient  feature  of  the  scattering,  the  Compton 
profile  ./(</). 

Simplified  forms  of  /.*  such  as 

/.*  > 2/.]  sin  \(j>  (52)| 

/.*  ~ 2(/.1/.2)‘  3 sin  \<P  (5  22) 

have  often  been  used  in  eq.  (5.19).  Use  of  cq.  (5.21 ) will  produce  an  error  in  q of  order 
[(/.2  — A i )./A j ] w hile  eq.  (5.22)  leads  to  an  error  in  q of  order  [(/.2  ) v.ij3.  Since  it 

introduces  a first  order  error,  one  must  be  careful  when  using  eq.  (5.21)  to  determine 
the  value  of  q,  especially  at  large  q values  and  for  the  case  of  large  scattering 
angles  and/or  high  incident  energy  photons.  Phillips  and  Weiss10  hate  calculated 
that  for  Mo  Ky.  radiation  scattering  at  90%  a q using  /.*  defined  by  eq.  (5.2 1 ) of 
095  corresponds  to  a q determined  by  eq.  (5.12)  or  (5.18)  of  10.  At  a q talue  of 
about  80,  the  difference  between  the  two  q values  rises  to  almost  10  per  cent. 
Williams11  has  pointed  out  that  the  use  of  eq.  (5.2 ] ) for  typical  X-ray  wavelengths 
leads  to  a scale  change  of  the  order  of  3 per  cent  and  hence  an  error  of  about 
3 per  cent  in  the  renormalized  value  of  J(0).  Front  this  discussion  it  should  be  clear 
that  such  approximations  are  not  suitable  for  accurate  work.  It  can  be  shown13  that 
q (eq.  (5.12))  determined  from  a non-relativistic  approximation  approach  at  the 
outset  differs  by  terms  of  second  order  in  (z.2  * ' \)t'\  from  the  q using  the  collect 
/.*  given  by  eq.  (5.18). 

In  conclusion  we  note  for  the  interested  reader  some  other  recent  derivations 
and  discussions  of  the  impulse  approximation. 13  15  Furthermore,  in  connection  with 
the  elegant  derivation  using  a time-dependent  operator  formalism  by  Fisenberger 
and  Platzman  (see  chapter  2).  it  has  been  pointed  out  by  Benesch  and  Smith1 5 that 
the  use  of  the  fundamental  assumption 

exp  ( - [/y0.  F]f3/2)  = 1 (5.23) 

in  the  many-eleetron  case  leads  to  an  additional  contribution  proportional  to 

2 dpi  dpi  F(k  f P|. p2  k|p,.p;)  (5(u>  - k2/l  k.p,)  (5.24) 

where  F(p,.  p2  j pi.  p/)  is  the  momentum  space  analogue  of  the  spin-free  two-particle 
density  matrix  in  position  space.  It  has  been  show  n by  Smith16  that  this  contribution 
is  negligible  under  normal  experimental  conditions  but  should  be  included  when  one 
integrates  d3<7  dQ  do  to  obtain  the  total  scattered  intensity. 


5.1.2  Accuracy  of  the  impulse  approximation 

Returning  to  the  basic  eq.  (5.2).  it  is  quite  clear  that  more  realistic  final  state  wave 
functions  than  plane  waves  can  be  used.  In  addition.  E j in  the  delta  function  can  be 
chosen  to  be  the  electron's  binding  energy.  With  these  changes,  analytic  results  have 
been  obtained  for  the  hydrogcnic  case.  Results  for  the  K shell  based  on  the  analytic 
expression  of  Gummel  and  Fax1 7 were  obtained  by  Fisenberger  and  Plat/ir.an6 
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q (a.u  ) 

F1GIRI-  5.1 

Comparison  of  M l and  IH  Compton  profiles  for  a K .he!!  electron.  7.  - 5.  7.  ; 70 

keV,  <$>  1 Sts  . Solid  line  Mi  and  da>hed  line  111 

and  studied  in  sonic  detail  by  Mendelsohn  and  Biggs.18  It  should  be  pointed  out  that 
Figure  3 in  Ref.  6 which  has  now  been  reproduced  in  two  review  articles1-4  is 
incorrect  as  it  shows  the  'correct'  results  lying  above  the  impulse  hydrogenic  results 
at  the  profile  centre  T he  error  was  fortunately  corrected  in  Fiscnberger's  later 
analyses  of  experiments  where  a Is  core  subtraction  was  needed.  A correct  com- 
parison'8 of  impulse  hydrogenic  and.  'correct'  hydrogenic  is  shown  in  Fig.  5.1  for  a 
case  similar  to  the  one  discussed  in  Ref.  6.  Note  in  this  fairly  typical  case  that  the 
curves  cross  twice.  We  shall  refer  to  the  correct  results  as  the  'Exact  Hydrogenic’ 
(I  II)  results.  Mendelsohn  and  Biggs'8  have  shown  analytically  that  for  Is  electrons 
the  impulse  hydrogenic  (II I)  profile  result  at  q 0 lies  above  the  FI  I result  as 

Ju i(0)  ./ i.tifO)  - 1 -t  0-l45/(fcoiio/Zf,)J  (5.25) 

where  k"  is  the  value  of  the  scattering  vector  at  q - 0 and  Zf,  is  the  effective  charge 
in  the  I state.  Cle  irly  the  correction  is  proportional  to  (binding  energy  i (momentum 
transfer)’  It  should  be  noted  that  this  statement  is  not  in  conflict  with  chapter  2 
wliere  it  is  stated  that  the  impulse  approximation  is  accurate  to  [binding  energy 
(momentum  transfer)  J:.  I he  latter  comment  was  made  with  respect  to  the  accuracy 
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FJGl  KE  5.2 

Compton  scattering  from  a 2s  electron  Z 5.  E,  20  kcV,  4>  ls0  Solid  line 
EH  and  dashed  line  JH. 


of  the  second  moment  of  St(k,(o)  while  we  are  discussina  the  accuracy  of  J(q ) at 
9 = 0. 

h'H  /.-shell  calculation  results  and  comparisons  with  experiment  have  been 
published  in  a series  of  papers  by  Mendelsohn,  Bloch,  and  Smith.19’21  These 
calculations  used  hydrogenic  bound  states  and  hydrogenie  continuum  stales  which 
are  represented  as  confluent  hypergeometric  functions.  Again  it  was  shown 
analytically 19  for  a 2s  state  that  the  IH  profile  lies  above  the  EH  profile  at  q — 0. 
Although  no  theorem  was  proven,  it  appears  from  calculations20  for  a filled  2 p 
sub-shell,  that  the  E.l  1 typically  lies  above  the  111  profile  result  at  q 0.  An  additional 
EH  /.-shell  result20  was  that  the  2s  state  can  exhibit  a secondary  maximum  on  the 
negative  q xide  of  the  profile  as  shown  in  Fig.  5.2  and  occurs  at  approximately  the 
value  of  q for  which  the  IH  result  exhibits  a plateau.  H is  secondary  maximum  is 
more  than  an  order  of  magnitude  down  from  the  primary  central  maximum.  It 
appears  to  be  related  to  the  structure  of  the  electron  density  in  the  original  bound 
state  wave  function.  Such  secondary  maxima  in  the  2s  slate  have  been  predicted 
recently  in  mu-atom  scattering  calculations.22 

Since  the  111  lies  above  the  1 II  i.-ult  at  q 0 for  the  2s  state  ..nd  the  IH 
lies  below  the  EH  result  at  q 0 fur  a 2 p state,  there  is  a cancelling  effect  in  the 
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calculation  of  a tola!  / -she!!  Compton  profile.  1 hus  the  F- ! ! correction  to  an  impulse 
to;  il  pt . is  than  the  corrections  for  the  individual  states.  I of  the  ca  -e  of  \!  o K i 
rad  i - ■ ■ ■ ; ugh  t70  from  tt  q 0 the  2s  IM  results  lie  6 0 pe  cent 

ah  ve  the  I tl  result--,  whereas  they  lie  7-1  per  eent  below  the  I li  results  !or  a 2p 
(avc  ug.-d)  election.  Since  there  are  6 2[>  electrons,  a 23V  per  cent  incrci-  • ol  111 
over  III  fo:  the  total  / -shell  profile  is  obtained,  lor  Mo  K radiation  scattering 
through  120  from  aluminium,  the  individual  2s  and  2 p corrections  arc  each  about 
20  per  cent  at  q 0,  but  the  total  / -shell  I'll  ,/i0i  lies  about  6 pci  cent  abo-  e the 
IK  value  In  section  5.6  we  will  compare  some  of  the'C  results  with  expcii  • .Us 
Clearly  /.-shell  corrections  to  impulse  can  be  quite  significant  at  convent;.'.'  •!  \ ray 
energies  and  scattering  angles.  Tabic  5.1  gives  typical  /.-shell  correction  - for  several 
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/.-shell  Compton  profile  maxima 
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Zfs(0) 
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JeM 

Jim  (0) 

J,hf(0I 

No 

5*99 

4-96 

2-609 

2- 54s 

1024 

Na 

6% 

605 

2*200 

2!  2s 

1034 

Mg 

7-96 

7 08 

1*908 

1-833 

1-041 

k\ 

8 91 

8 II 

1*692 

1615 

1-04S 

Si 

9 88 

9 11 

1*521 

1444 

1053 

Ar 

13-69 

13  08 

1077 

1020 

1056 

elements  for  Mo  K y radiation.  As  one  might  expect,  the  corrections  to  impulse 
increase  with  increased  binding  (i.e.  larger  atomic  number). 

Since  the  EH  approach  to  the  correction  of  impulse  profiles  utilizes  a screened 
hydrogenie  model  for  both  the  ground  and  continuum  states,  both  electron  exchange 
and  correlation  are  neglected.  The  use  of  a more  realistic  model  for  the  evaluation 
of  eq  ( 5.2)  has  been  considered  by  DeCicco  and  collaborators’3"4,  and  Grossman 
and  Mendelsohn.’’  These  authors  have  used  an  atomic  central  field  model  potential 
due  to  Herman  and  Skillman  (H$).’f’  It  employs  the  full  Coulomb  potential  with  the 
Slater  Kvc a!  approximation  to  the  non-local  Hartree  Eock  exchange  and  the  so- 
called  Latter  r ' 1 tail.  It  corresponds  to  the  choice  of  x 1 for  the  exchange  parameter 
in  the  Hartree  Eock  Slater  ( .Yy)  model  approach  of  Slater  and  co- workers  vv  hich  will 
be  discussed  in  section  5 3 for  impulse  calculations. 

In  their  method  for  calculation  ol  exact  Herman  Skillman  (El  IS)  Compton 
profiles.  Grossman  and  Mendelsohn  use  a modified  version  of  the  HS  continuum 
wave  function  code  of  Munson.’  1 hey  have  found  that  of  the  order  of  2u  or  more 
partial  waves  must  be  included  to  obtain  convergence  to  three  significant  figures  for 
a typical  point  on  the  profile.  The  method  has  been  shown  to  give  agic  meat  with 
ti  dytic  EH  results  to  more  than  three  figures  when  a hydt  g nic  fathei  than  a HS 
potential  : . used  This  method  can  be  used  to  calculate  the  Compton  profiles  ol 
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electrons  mi  an\  '■  cl  I !’•  ■'!>  I I IS  t .V  alations  were  c.u  lied  out  fur  the  iu.n  ease 

: , . . rr.ii 'noil  ,ibuu  and  wci  c h 1 ml  U> . i g i .1 ! i .•-!  with  the  cxp:-;  i mental  r e-i  Its  o > 1 1 the 
whole  piolile.  In  addition  an  • ; :Ne  i ah  elation  w as  pcrfoimcd  using  a e-'  • ^ 

stale  1 IS  w.no  funs' m i.  \.  . in  I '..  II  or  reel  huts  to  impulse  of  ah< 'tit  3 pot  rent  wore 
i>h;..ins-d  at  ij  l).  I Inis  the  i n . i . ■ li'nle  of  the  corics  lion  fur  neon  seems  relatively 
iiuleji.  idem  of  the  Iheoretioal  bound  Mate  wave  function  approximation  used. 
\nothei  iitlei csting  preliminary  re  alt  . 'htaine.l  w ith  the  partial  wave  IIS  method  is 
ill, it  fur  60  keV  photons  s.  atl  dug  from  krypton,  the  corrections  to  the  IHS  ./(0| 
nil  aes  for  the  indi\  idu.d  orbitals  may  he  of  the  or  der  of  a few  percent.  I he  cun  act  ion 
to  the  impulse  J(0|  result  for  the  entile  atom  is  significantly  less  than  this  hei.iase 
of  s ,i ui citation  edicts. 

\s  discuss  -d  in  the  next  section.  : npulse  I lai  tree  l ock  Slater  (,\y)  j utiles  are 
si:  si.ne  to  the  ehuice  of  the  e\cl  mge  parameter  z.  especially  for  outer  orbitals.  1‘he 

nc  ; _ it  of  the  I IIS  neon  calculations  with  expsaiuicnt  may  l.tdieatc  that  the 

’ ■ re  of  j.  1 iii  ; v 1 a\e  sp  . ' il  advantages.  Such  a eluie.j  of  y is  known  tii  j 11  in 

tli.  j )••.,!  w.ne  function  loo  much  as  compared  to  a Hartrce  lock  wave  functiun. 
In  the  l!!S  c.de  il.  liun.  this  v.a\e  function  chi  ice  may  con  j ensile  in  j art  for  using 
the  rite  potcrdial  for ealeuhiting  the  hound  and  continuum  wave  functions,  ti  e.,  the 
fi  o/cn  1 1 re H u - i'  iiion |. : 'd  in  part  for  m oi  ing  con  elation  and  rclali'  istic  effects 
\ di  »,_■  nssiuti  of  the  ‘fi  • ■ n io  e a ppi  oximat  inn’  has  been  given  by  Cm  rat  el  al.~ ' and 
h\  Smith. A further  possible  '.mpru\i : lent  ter  the  I IIS  method  could  he  the  use  of 
a more  realistic  value  of  y.  such  as  yt  P deti  r mined  by  Sabin  ml  Si  null  k and  disc  assod 
in  section  5.3. 

Il  has  been  asserted'' that  the  impulse  appro-in  atie  -.gives  very  a«.c  mate  le'ults 
fur  large  q since  the  electron  is  then  ejected  w ith  high  mo.mntum  . d the  p’.me  wave 
representation  becomes  valid,  this  wi  aid  certainly  be  the  case  if  this  wee  the  < tly 
approximation  inherent  in  the  derivation  of  the  impulse  a p,  loximation.  In  the  delta 
function  the  hound  slate  electron  lias  also  been  taken  as  fiee  and  it  is  not  clear  lic'W 
this  affects  the  asymptotic  character  of  the  results  far  out  on  the  profile  Figure  5.1 
for  a K electron  shows  that  the  FH  pr  >,:’e  lies  siih'ta.:  dally  above  the  111  protile 
if  one  goes  out  far  enough  in  q. 


5.2  MOMKVri  M DISTRIBI  I IONS  \ND  I Mi’ l I SK  ( OMI’ION 
PROFII.KS 

In  the  previous  • tion.  we  have  dis  ssed  the  cot  lection  between  the  impulse 
Compton  pi  file  1 ction  J{q),  and  the  io  icnluni  distribution  within  the  context 
of  a one-electron  model.  As  shown  by  Bcacsch  ..nd  Smith.1'"’11  ‘ one  may 
generalize  the  relations  (5.14)  and  (515)  to  reflect  that  J (p.-l  and  ./ (</ ) arc  one- 
! tsional  p < ieetions  of  p(p)  p(p|  p)  w here  p(p  | p’)  is  the  one  par  tide  ch  rge  (or 
'isin-frce)de’is:tv  matrix  in  itk  pent  uni  space  for  the  \ -elect: op.  sy  stem.  It  is  d 1 d 

ill  : u i.  r.  sj  :,ce  / ( A •.  \ . . \ \) 

momentum  and  spin  coordinates  of 

electron  j. 


by  eq.  (5  '61  in  te;  is  of  the  wave  function 
where  \\  ip^.  s_,l  denotes  the  combined 
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;( i»  ‘ p'>  I **(A..Y. V,)/!'  V. \\l(d.Y).  (5.26) 

It.  tl  . |d.Yh  denotes  the  combii  l operations  of  intej over  the 

i.  i coo  linat  s and  su  tion  over  the  spin  variables  of  electrons 

2 5.  . . \ Jii't  as  the  wave  function  in  momentum  space  is  related  to  the  wave 

funs’’ m in  position  space  via  the  Di  ae  i cn:rier  transformation,  so  is  yip  p’l  related 
to  its  po  iti  at  s,  ; • analogue  ;•( r r")  by  *' 

yip  [ p ) - (?*)  3 j y(r  | r')  e •>'  r dr  dr’  15.27) 

Thus  the  fundamental  relations  are 


J{pz) 


PlPx-PyP;)  dPx  d/’j- 


JUj)  : P 'l(p)  dp 

J i« 

where  the  radial  momentum  distribution  /(/’)  is  defined  bv 


/(/ ) p- 


pip)  df] > 


(5.2-S) 


(5.29) 


(5.50) 


These  relationships  clearly  indicate  that  knowledge  of ./(;/)  leads  to  info:  nn.  >tion 
ahoi.it  the  momentum  distribution  However,  it  is  logical  to  ask  at  this  point  why  we 
should  he  interested  in  the  momentum  distribution  except  that  it  is  accessible  through 
the  impulse  approximation  from  the  experimental  Compton  line  shape.  Can  it  add 
to  the  knowledge  of  electronic  structure  which  we  can  derive  from  the  familiar 
chaige  density  p(r)  - y(rjr).  itself  derivable  from  the  experimentally  niensu.-ouMe 
atomic  fori  i factor  /’( s) 


/•Is) 


p(r)  e" ' ’ dr 


(5.5 1 ) 


‘I  his  question  is  best  answered  by  considering  Fig  5.5.  where  theoretical 
c:d  'i’at’n ms-"  of  /(/>'.  J(i/).  / |s).  and  D[r)  arc  presented  for  the  1 h state  of  atomic 
va.  ad.-.r  1 1 s ’2' '2/«''5s '5/ ''5i/'4.. ’).  Here  D\r)  is  the  radv!  ch.uge  density  defined  by 

D{r)  r2  j f>( r)  dC>  (5.52) 

1 1 compl  ;nt  y and  invet  nature  of  the  inf  mat  n presented  by  JUj)  and 
/I/’t  a«  1 by  /’|s)  and  D\  ) i cl  irly  seen  It.  the  plot  of  D[r)  we  cai  ol  c the  pe  tl  s 
cor  cspouding  to  the  K.l.  and  \/  shells  bet  lose  sight  of  the  V shell  in  the  tails 
of  the  other  shell  Ot  the  othe  hand,  the  radial  mom  turn  density  Up)  curve 


p ( a u I 


c/  (a  u . I 


ic) 


< d I 


\ IGl  RH  5.3 

of  Harircw  F.vk  ;«.vd  *.  .!  ' • ‘ :i"  for  •.  ' ,..n  \~F\ 

(a)  Radial  c)  gc  density.  i)(r).  (b)  Atomic  form  factor,  fish  (cl  R ' i 

• 5 /:/ 1 »d)  T npuNe  Coir.pt©:.  podi u*.  ./(</) 


‘•'iou>  the  V ;nid  A/  shells  res-  ccti\ ely  . -s  p increases  while  the  I.  shell  is  se^n  as  1 lie 
shoulder  on  the  cane  at  6 ;ui.  and  tbc  K shell  is  lost  in  tlic  tails  of  die  others. 
C ■ 1 f the  _ tail)  ble  q ntities  lit  les  that  1 

1 dir  tlv  ’ 1 : e while  the  t • 1 1 ’ es  in  tl  lope  of 

Jhj)  d,cn v.  tsh  the  \ and  \J  shells  uhtlc  the  core  d< '-ninales  the  tads,  f he  >Vpe 
o f Jhj)  is  sirnplx 

J'Ul)  yj  ’/(,/).  cv33) 

It  is  dear  tl  .it  the  eh:-’:'e  am!  • • s;.vim:!r>  d; nsilies  • !*»  j o\  ide  <.  v , \ t :r\ 
r iition.  One  of  tl  ! : s of  the  C plon  s<  ' ique  lies 

fuel  i • l 1 he  i . ' ■ • ; i ’ e ' 1 t La  ' ■ 1 1 ^ is  ' ■ a * e,  f-  ■ ■ : t ] # ^ i . y i . . . l 

the  cv  j l ; e>. ' . t - .i : e '<•  'St  neei ; i .ii e ( ala  l h.i pier  3 1. 


aiom>  in 


Anotl er  .I-,  r of  It  s i el  ti  >n  betw  ■ F(s)  and  ./ ( j)  w is  shown  by  Be  ;scl 
alt.  and  Si  : ;t '' , ' * n.'  ')  tint 


Thus  F\ s)  is  a convolution  of  the  entire  one-partiele  density  matrix  in  momentum 
space  ylp'p-)  whi'e  ./(./)  is  a one-di  nensi .anal  projection  of  its  diagonal  component 
ftp  I P)-  The  stronger  smoothi  g nature  of  the  co  ti  > in  ivexpl  in  why  correla 
tii'n  correction » to  atom'  fort  a factors  arc  UmI  illy  si:  tiler  than  those  for  Couq>s..u 
profiles  which  in  turn  are  smaller  than  those  for  th.  momentum  d: -tribuii  ni 
Inspection  of  f ig.  5. a cot. firms  this.  We  shall  return  to  the  subject  of  correlatin' 
corrections  in  section  5.4. 

It  should  be  rec  died  that  if  one  proceeds  from  th.  position  sp  • .e  rej  . cs,.  ,t  ition. 
deternii  .ation  of  p(r)  and  F( s)  requires  only  the  ding  u Component  ;ir  rl.  whi’e 
eq  (5.2  7)  shows  that  the  full  -;(rj  r‘)  is  required  for  j(p|p). 


5.2.1  Calculation  of  p(p)  and  ./ ( c/ ) for  atomic  systems 

The  preceding  discussion  makes  clear  that  one  way  of  calculating  pip)  and  ./Ip) 

i>  to  use  the  prescription  (5.26)  to  calculate  ;((>  p)  from  /(.Y,..Ys A\).  In  turn. 

the  ware  function  in  momentum  space  may  bool  '.lined  by  I brae  Tout  ier  ti  a reforma- 
tion of  the  position  space  wave  function,  i.e  . 


/«  Yj . A\ A\t  (2a) 


t/'(.v1..Vj V.v)  exp 


- < i P;  ■ r>  ) 

j - l / 


where  \y  (rr  st)  is  the  combined  position  and  spin  coordinate  of  electron  j and 
(d\)i  has  the  analogous  meaning  to  (d  Y)>  defined  abuse.  Since  many  position  space 
wave  functions  arc  available  fot  atomic  systems  or  arc  easily  calculable  today,  this 
transformation  is  the  commonly  followed  procedure  (or  the  equivalent  transforma- 
tion (5. 27)  of /(rV) ) and  will  bedisais>ed  in  some  detail  later  in  this  section  How- 
ever by  writi.  g the  Schrbdinger  equation  as  an  equation  in  terms  of  /.  it  is  possible 
to  determine  / directly. 

I . Dover  so/ti.v'i ’ii') J]ir  /.  1 wo  approache - may  be  conk  ' med  1 he  first  is  a differential 
e ; ■ . i in  momentum  -j  ce  where  a simple  replace  lent  of  x by  id  dpx.  etc . is 
nude.  Since  the  r ' potentials  involve  differential  ope  itor.s,  this  approach  is  not  a 
simple  on:  Hyllct as-.'''  did  solve  this  differential  equation  for  / in  the  case  of  the 
hydros  i atom.  In  the  second  approach  application  of  the  Di  ic  Fourier  transf 
tioti  to  the  position-space  Sehrodi:  _•  -r  cqnatio  • lea.. Is  to  an  it. ’eg  a!  equaimn  for 
A-  As  an  exan  pie,  it  takes  the  following  fot  i fot  the  one  electron  hydrogen  ic  atom  : 


2/l/ip) 


p p ’/(p’)dp'. 


t'a  ill;-  case  it  .v.i,  solved  exactly  by  Took.'  " Hi.  analogous  equation  for  the 
twi  o . problem  vv.i , at'aeked  by  McWeeny  and  Coulsi  n’  ’ and  by 


Hculeis. >n  and  's  -hoi r"  ’ • ng  ,<n  itc:  dive  method  of  vilti|ji-ii  I 'c  chi’iceltie  • ' ch 
thev  encountered  ..re  • •;ic  t for  the  n.  gUvl  (’fils  i;>c  ■ v then. 

2 ' ' Fourier  in  union  f t/r  or  ■ Although  the  di  reel  sol ui ion  n < 1 1 1 > 
im.  ill  today  due  10  the  > -..plcxities  involved  in  •n..t.ing  it  cotnpur  ti< • 'illy 
t r.i hie.  \1c\ki.cny  and  ( 1 dson  y w ere  led  to  o 'd.-r  it  ,0  a 4 ^ ; 'c 

alternative  to  the  Dirac  Fourier  trnn.sfo* i.'.ation.  fliey  c«>  '.’’Killed  that  ‘...,t  j > •.  as 
impossible  to  transform  to  momei  turn  'p.ve  the  most  arcuate  wave  fin  c:>o;:s  in 
which  electron  eorre’ntion  is  . ’equately  represented.  the  inlejials  involved  K ig 
quite  ittsiipei ableV'u  Ind.-ed.  be- :de  the  difficulty  of  the  integrals,  it  is  cei'amiy  a 
l.iborious  undertaking  to  ke  the  total  wave  function  in  position  space,  coi  tpete 
its  3.Y-di  uenOua!  Oiiae-I  ier  transform  (eq.  (5.35) ).  squme  the  tr..-.sfoim.  ; d 
then  integrate  out  the  momentum  space  coordinates  of  the  .V  1 elections 
(eq.  (5. 16) ).  In  the  case  of  an  independent  particle  model  wave  function,  i.c.  one 
which  is  icoreseni.'.ble  bv  a 'in  tie  5!  -er  deter ' inant, 

tjij  x 1 ) t fijxj)  ...  >!'J  v.v  I 

: l/'h(V'l)  (,'/•.(  Vs)  •••  (M-Vv) 


'!>,(*  i)  (,'/,( Vs)  ...  ^,(\\) 

one  needs  only  to  transform  each  individual  spin  orbital  by  means  of  eq.  (5.‘>|. 
Similarly  in  the  ease  of  a wave  function  which  is  evpanded  in  terms  of  a set  of  spin 
orbitals,  the  transform  of  each  individual  spin  orbital  leads  to  / expressed  as  the 
s:  me  exp;. nsion  of  Slater  dele  u iiiar.ts  bat  oxer  the  trarc-foi  .vd  set. 

An  alternative  appn  ch  de\  loj  ed  by  13  esch  and  Smith  5 allows  . turn 
distrib  a tioiis  to  be  obtained  fi.mi  the  r..ost  on  ■ ’ic;.u  d posit '•  <n  s;  ■■  ce  wave  r..:  >c:i  ms 
wherein  elect ron  correlation  is  adequately  represt  ited  Since  the  < le-particle  cl  ge 
density  matrix  */f  r j r’i  may  be  expanded  in  terms  of  its  c iger  functions/1-  ' the  natural 
orbitals  (NO's)  </q{ r). 

;(r]r’)  - (,y,(rti,'/*(r'l  (5.3S) 

v here  is  the  occ  uj  ■ tic  n number  of  yqlri.  k iowh  dge  of , (r  ’ r'i  in  tci  ms  of  its  NO's 
enables  one  to  im  eq.  (5.26)  to  obtain  ftp'  p‘)  exp.  tded  in  ' ns  of  the  motncntuin 
space  NO's. 

Vfp ! p'l  - V/./  i'l/Np'l.  (-N39) 


y'(Yi.v2 vj  (.V!) 


Thus  in  the  case  of  a ytr'r'l  of  finite  o ie-i.mk  (i.c.  th.>  ’ er  of  nonzero  in 
eq.  (5.38)  is  finite)  one  ne  ed  only  o ipnte  the  13"  ,<c  ! ■ , Tei  t...  -f.  m (eq.  (5.9) ) of 
each  individual  NO.  lids  method  has  beer.  ;•'.••!  to  * ’ ’ V'  .'  ' ip  TVs  by  a 
number  of  null  nrs. 

Hie  c.l'C  of  ave  f'.uvt  .■  .ns  ; 1 ( \ 1 • ( * . . ' ■ . , iiv'.e'ler- 

electronic  coordinate  rt  is  interesting  fhcic  ..re  two  ipj  > -ches  w!  eh.  ay  be  u-  d. 
One  can  ....  nd  the  n.  t etuis  in  a set  of  « e : article  f.  net'.  ' w’u.h  in  o.avticc 
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are  a truncate  ! set  43  The  NO  expansion  is  then  available  and  th  !:.  . ! Stuit! 

procedure  is  applicable.  Alternately,  one  may  pro»  ’ t tal  late  y(p  j p)  d ■ a s 
shov  n re  ■ tly  by  nhui  " an  i Ben  h ' Sot  eoft  Is  involved  a 

si  ii pi  i and  require  extei  sive  eoi  puter  ti  e is  >h  wvj  n tl  cal  la  of  B I 
for  Ho-1'  I ..  o earlier  calculation'4  ■•4  of  the  ti  • me,  • a d;  ' !ro 

Hvlleraas-type  function'  for  atomic  helium  ha v e bcei.  rcp>»  •,  .1  Ore  of  ti  . •4*' 
contained  errors15  while  Hicks*7  chose  wave  funetu  is  wl  h incfud  i only  terms 
of  the  form  rf,  whereas  rl2  terms  are  needed  to  [ ly  re|  t sl< 

elect 'or  cusp,  llie  use  of  r\,  tei:  >>  made  tl-  ex  ’■  of  t'  n 

distribution  much  simplci. 

3.  Dime  Fourier  tmusf  >rm  Uion  of  orhiu;!*.  With  the  e'ce;  • • . ot  it  ea-e  ol  wave 
functions  which  explicitly  contain  the  intcrc'octronie  coordit  a'.  % riJk  we  have  seen 
thattheDirac  Fourier  transformation  of  the  tV-electron  w a eft  :ti<  i^ortheone- 
p ; •!  • density  matrix  ; f r ] r ' ) ultii  itely  rests  on  the  ti  ol  orbital 

flier- 'i ins  from  po iition  to  i:»  tin--  spa e-a  fee  1 5 9)  ).  In  the  s I tp|.. a o.i-onfj'  • ic 
hyd  ogen.  Pod  >lsky  and  Pauling43  transfoi  d the  positii  space solut 
ol  bo  tie  exact  bound-state  y(p)  with  which,  the  later  direct  e:i!-:til  itionx  ot 
Hyllera  is-*®  and  Fock  ’ wee  ir  agreement.  The  solutions  which  cot  respond  to  the 
fainih  ><■  po'!*;..  , spare  solutions 

i j/(r)  ■ J„;(r)Ylm(0,<p ) ( 5.40 ) 

are 

-/_{[>)  ujp)  Y:J0P.  <!>,)  <5  41) 

where  <p)  and  Y,J0  tft  ) are  the  usual  spherical  harmonics  in  position  and 

ii’.oa  ei.fii  i space,  respectively.  This  illustrates  the  isomorphic  nature  of  the  Dirac 
Four'.,  trai  sfo  mution  (eq.  (5.9)).  namely  that  s-otbitals  transform  to  s-orbitals, 
p's  to  p s,  etc.  Instead  of  writing  down  the  general  solution18-*'1  which  involves 
Ci.’g. ■ poly nomials,  we  list  the  Is  and  2p  solutions  for  future  discus*;. -n  and 
refe  race. 
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(32/ 
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w’- . , / I*  the  tiuclc-tr  ch  mge  and  / 7.  2. 

In  ti. .-  c .sc  of  i many  e!r  -<ron  system,  we  assume  that  an  orbital  y/r)  such  as 
■'  n il  : ' orbital  or  l|..r*tee  I oc'\  orbit. d may  be  wnUe.t  m the  position  space 
i eq  1 5 40 1 or  in  the  r , > r a gene:  .d  case  as  a linear  con-Si;  ition  ol  such 
re  -tr'ete  !'  fm  ..mins  f or  such  atomic  orbitals  the  radei!  funct  ' n /„;(.••)  is  U'.:  b!y 


ws.tt.-n  as  .:n  <.  . ton  over  i rime  set  >>'  . • • y l ic  i -ions  or  ;■>  a i 1 it 
tn’".d  iiioa.  rite  i>i-  .'nlum  spa.  e othital  i-  of  the  for m of  eq  (5.41 1.  ; e..  / jh  ' .ts 
the  ' • far  dep.  d ce  as  ^ r)  and  the  t •’  1 f u • • • ii  . i by  lie 

v\  i .v-tr>  -di p dent  tr:  .'sfo'  • ".aiion. 

'■■■  Ip)  I i)\  2 r.  I ],[,,) f. iir 

• 0 is  48! 

( i)‘p  '//,,(/.) 

•.»  I.cit.'  I he  I fardel  Ira  i sforin  f /„,(/')  involves  a q ’ . it  a)  H.  e)  function  j,(r)  of  the 
order  of  the  en.idar  moment  i in  quantum  number  / In  this  case  it  the  H;.  Vcl 
, nsforiii  of . dcr  / J,  1 of  the  f.  nction  rf,Jr),  i e., 

//„.(/>)  I [r/;,(r|]./,  j ;(/’.  )(/  / r ’ dr  (5.  ','?l 

Jo 

where  J,.  i 2 is  the  ordi  .ary  Bessel  function.  For  a mi;  tl  .i  of  p He  ft  .ms  for 
the  u„i(pi  hate  been  derived  in  the  literature  for  various  sjv.  .•!  and  arc 

discussed  and  tabulated  q ;ite  generally,  U-icther  with  the  , d ./(</)  In  S.  . t h 

and  K ii’>er.';“  For  atomie  ealeulations.  the  o ...  ly  '.>ed  f.  ' d ft  ns  for 

/n,(r)  are  the  hvdiogi  r .ic  f.. notions.'*'  the  nodclcss  Sitter  type  oil  '..Is  (S  I ()) 

fjr)  [Cz)-"- \lC,r  Jij'-V  ’e  " (5o0| 

and  the  ('.:m  inn-type  orbitals  (<  i I O) 

./„,(/•)  - [.ICaf  1 ; F(n  i)J*  V ' e (5  51) 

In  the  ease  tif  orbitals  which  a e either  tabulated  nt:  i ic  dly  or  for  which  the 
analytical  transform  is  not  no. liable.  ■ a nciic.d  t.-.l  ,i>es  t 1st  be  o,  ; toyed.  I he 
accurate  numerical  evaluation  of  Fourier  sine  and  cosi  c Irai  • - h rswel  ive 
in  eq.  (5.4S).  is  beset  with  dillk ulties  due  to  the  . pid  . o ''  lion  of  the  i t ma  ,d 
As  p - r. . ,-ine  and  or  cosine  functions  oscillate  ...  1th  such  • idity  that  o n in  the 
ease  of  radial  functions  f„i(r)  that  very  rapidly  ippro  h Atro  as  r apj  ' $ 
infinity,  the  transform  is  made  up  of  an  evtr.incly  lame  number  of  'J  . dkant 
positive  and  negative  comributions  of  nt.-.rly  equal  si,c  'hat  iisc  ii.rn  the  cyek*s 
of  the  tri  nometric  functions  involved.  It  is  for  this  son  that  the  rt  Its  pri  ted 
in  ceral  cent  st  dies  involving  nu  ericn  1 1 fi  ti»  of  ati  ii  bital  ■ e 

hr. h inc  in  the de'ii ed  at'euau y flteaect  ••.•locvah  vf  these  t>  sf< ■.  ms  1 as  1 n 

considered  recently  by  11  d t.ar  and  S • 1 1 h.' f’  by  B-.  i-eseh**  a id  by  P dsirup‘:  ■ th 
tl  e f<  • -i  : 'tl  ors  mpl  tsi  ig  tl  . . valuati  n to  ...  • in  lied  acy. 

4.  I.. t < ; /■.  . : ■ ‘ ■ s , ,r  / / .1  ( ; | ( ' ; in  i cuts  of  the  m 11 
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a , ' p|  . i|  ■ _,t  Tl  best  kn<  is  p'  which  is  twic  th  ki  tic  energy 

bv  vii  ; of  tl  • I tl  s just  twice  the  negtriv  >f  the  total  at 

I , ' • : ! deliniti  eq.  (5.2‘i),  shows  that  1/2  l2  is  ju st 

i tl  peak  val  Compton  profil  wl  ile  pn  is  .V  t 1 o{  it  p* 

ar' ...  i . t!>  • *o  .■>;  > I'll-  re!  >’  .isve  Hiet  equation  written  in  the  Pauli  apptoxima 
tie>  ‘ 1 |t  • iti  t to  : that  the  possible  integral  moments  p*>  are  confii  1 

t..  th.  r.i  2 < a 4 fot  the  exact  wave  function  am!  unit  tot  0 < /t  4. 

:ie  lKc/">  (5v*) 

u'  . ,/•>  i>  the  nth  mot  ,e  t of  the  Comp:  >n  profile  function  ./(«/». 

V>f  | q"J{q)  d,/. 

lo  u Jet  stand  these  st.  ?en;ent»  it  is  neeosx  ->  to  invest^  the  sma!’  ard 
t*  • .•  p be’  . - of  /(/■•)  |or  ! pi ) and  si  mil  • R the  small  and  'em-  q behavt*  ••  r 

of  }(./) 

1 he  expansion  of  Jiq)  about  q 0 is 

jit;)  j(0)  Ttpiou;-  i^v'iojt/4  •••  ,5  ;S| 

a'  sin  " rly 

I(p)  4r;p(0)/h  4 j,-p"(0)pJ  f (s 'ol 

Si  ; thee  vpet  imental  measurements  are  the  most  accu  te  around  the  profile  pe  4 . 
B,.;  h .-ii'"  sue  wee.'  that  />|0)  and  /V|0)  can  be  extracted  from  experimental  ./(»/) 
c ves  and  compared  with  theoretical  calculat  is  Use  of  eq  (5.56)  in  eq  (5.52) 
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5.3.1  ! ' r ■ r . c Koek 


Since  it  is  ot  possible  to  solve  the  Schrddingcr  c < i exactly  for  systems  e<  in- 
ing  more  than  a single  electron,  approxii  te  rhods  for  the  d<  ition  of  the 
wave  function  are  used  which  i e is  ally  b sed  o the  < erg)  variati*  il  pri  > iple. 
1 fe  si.  ip'cst  • .at/  :s  that  of  the  s elc-Sh/er  o’  4ei  ..nt  defined  in  cq.  (5.37)  or  a 
li.-c«ir  combination  of  such  d-.Ni  as  chosen  so  that  the  total  wave  function 
poss  ■ ses  the  proper  symmetry  , roperties  (L2  and  S2)  of  the  state  in  juestion.  The 
i '•  ita’s  b..olved  a ill  be  s.  'hit ions  to  the  at  m ic  1 1 K ..dons 


jV2  £ r ■■ ■ i;  I, 


l5fO| 


where  l'c  is  the  Coulomb  pole4 it  il  due  to  the  e!«ct ionic  charge  /<(r).  is  the  non- 
1 ■ d c' change  potential,  arid  the  orbitals  </• , ; r ) . a . . nod  to  be  C't:  ict  d to  the 

4 1 .'.  10)  Since  Hc  dept  ids  on  the  N o upied  orbitals  [^(r)]  ’ solution  of  tl  ,-se 
...  4 ns  was  foi  lanv  years  a stumbling  block  in  the  tl  ! 4 ■ ipti  .of 

.4o  ic  st nict ure.  As  .-.  ro  It.  tin  ear  . vi  v .'  :v  4 r"  of  (’  ( ■ n J . Vs 

for  ii  .::i) -c'c-. iron  atoms  cl  .ployed  as  unde  . , , . ■ :io  -s  to  the  - ’ tions  of 
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FIGl.'RK  5 4 

Radial  momentum  density  / ( p) / .V  for  the  ground  states  of  the  first  row  atoms 


eq.  (5.62),  either  hydrogenic  or  Slater-type  orbitals  with  shielding  factors  chosen  on 
the  basis  of  Slater's  rules68  or  similar  criteria.  However,  many  of  these  calculations 
were  reported  with  7(0)  fixed  to  a common  value  rather  than  the  profile  area 
normalized  to  (V/2  as  expressed  by  eq  (5.52)  for  n = 0.  The  first  extensive  study 
was  made  by  Duncanson  and  Coulson6*1'70  who  reported  calculations  for  the  neutral 
atoms  from  hydrogen  to  potassium  (Z  = 19)  using  STO’s  (eq.  (5.49) ) or  with  orbital 
basis  functions71  of  the  Morse  Young  Haurwitz  type.  They  were  the  first  to  observe 
the  momentum  space  shell  structure  illustrated  in  Fig.  5.4  for  l(p)/N  for  the  first 
row  atoms  (Li  Ne)  and  which  we  discussed  in  connection  with  Fig  5.3  for  V.  They 
noted  that  s-electrons  give  narrower  or  more  compact  / (p)  distributions  than  do 
p-electrons  and  were  thus  able  to  explain  via  the  aujhau  principle  the  successive 
expansion  and  broadening  of  the  Up)  curves  observed  in  Fig.  5.3. 


Table  5.2  Extensive  tabulations  of  HIE  Compton  profiles  for  neutral  atoms 


Authors 

Wave  function 

Atoms  1/1 

Relativistic 

Weiss,  Harvey, 

Clementi72 

Li(3)  Ge(32) 

No 

and  Phillips'4 
Benesch*0 

(STO  basis) 
Numerical 

As(33|  Ybi 7<)| 

No 

Bums.  Mendelsohn. 

Numerical 

Hill  Kr(V>i 

No 

and  Mann 

Biyjs.  Mendelsohn. 

Nu  menc.il 

Kri'N  NbtlO?) 

Yes 

and  Mann"-" 

Biges.  Mendelsohn. 

Numerical 

Hel2l  Nbl  102) 

Yes 

and  Mann'”’ 

With  i he  advent  of  modern  computers,  it  h;is  become  possible  to  solve  eq.  |5.(>2l 
by  the  self-consistent-field  I St  F)  procedure  either  in  terms  of  analytic  basis 
functions  * or  numerically,  ' As  a result,  impulse  Compton  profiles  of  Hartree  luck 
quality  have  been  calculated  in  recent  years  by  a number  of  authors.'0  81 

The  most  extensive  and  useful  tabulations  M 80  are  listed  and  described  in  Table 
5.2.  The  Compton  profiles  for  the  numerical  Hartree  Fock  wave  functions  should 
be  considered  as  a reference  standard  for  comparison  of  calculations  made  with 
either  more  elaborate  wave  functions  or  basis  sets  intended  for  use  in  molecular 
calculations.  In  the  first  four  columns  of  Table  5.3  we  have  listed  J (q)  for  four 
Hartree  Fock  quality  wave  functions1'-81,82  calculated  by  Smith,  Brown,  and 
Beneseh81,83  in  order  to  indicate  the  effect  of  basis  set  quality. 


Table  5.3  Comparisons  of  \arious  correlated  and  Hartree  Fock  impulse  Compton  profiles  for  boron  ( /')' 


Compton  profiles  were  calculated  b\  Smith.  Brown,  and  Benesch. 81 
dementi  " ‘ 

Bagus-Gilbert 1,2 

Numerical.*1 

2s-2p  Correlation.1 1 1 

The  densrfs  matrices  of  (he  Schaeffer  Harris  functions' 14  arc  from  Brown  and  Smith 


5.3.2  Local  density  models 

Because  of  the  complexity  of  solving  the  Hartree  Fock  equations  (eq.  (5.62) ) which 
is  especially  due  to  the  presence  of  the  non-local  exchange  term  Vx,  there  has  been 
a long-standing  interest  in  the  development  of  local  density  models  for  the  description 
of  the  electronic  structure  of  atoms,  molecules,  and  solids.  Because  the  atomic 
Hartree  Fock  equations  are  solvable  as  a routine  matter  today,73  atoms  represent 
an  ideal  lest  case  for  such  models. 

The  Thomas  Fermi  (TF)  and  Thomas  Fermi  Dirac  (TFD)  models  have  been 
applied  with  some  degree  of  success,  especially  for  heavier  atoms.  The  corresponding 
momentum  distributions  and  Compton  profiles  have  been  considered  by  several 
authors. 05,84 ' 80  The  TF  model  gives  a profile  of  the  wrong  shape.  The  peak  value 
is  quite  high  in  comparison  with  HF  even  for  a heavier  atom  such  as  krypton.  The 
value  of  J( 0)  may  be  shown  to  he  CZ1  3 where  C is  approximately  equal  to  3-2. 
This  wrong  behaviour  for  small  p is  due  to  the  r 0 asymptotic  dependence  of 
p(r)  in  the  TF  model  instead  of  the  proper  exponential  decay.  On  the  other  hand 
the  TFD  method  cuts  p(r)  off  after  a finite  distance  and  hence  leads  to  a more 
realistic  behaviour  for  small  p and  q. 


mm. 


q la  u.) 

HKb 

HFb 

HFC 

HF“ 

2C.T 

Pol.' 

FO' 

I87CI* 

00 

29924 

29887 

2-9889 

2-9896 

29507 

29685 

2-8878 

2 8954 

01 

2 9135 

2 9107 

2-9110 

2 9113 

2-8757 

2 8911 

28187 

2 8250 

02 

2 6918 

2-6911 

26912 

2 6912 

2-6648 

2-6735 

2-6224 

26256 

03 

2 3684 

2-3697 

2 3695 

23692 

2-3545 

2-3561 

2-3314 

2 3310 

04 

I 9984 

20006 

20002 

20000 

1-9964 

1 9929 

1 9918 

1 9885 

05 

13118 

1-3129 

1-3127 

1-3128 

1-3236 

1-3173 

1-3413 

1-3363 

08 

08432 

08429 

0-8432 

08432 

0-8574 

0-8527 

0-8779 

08750 

to 

05788 

05784 

0-5786 

0-5785 

05903 

05873 

06056 

06050 

As  an  outgrowth  of  the  ideas  of  the  Thomas  Fermi  model,8  there  has  been  a 
great  deal  of  attention  paid  in  recent  years  to  developing  methods  wherein  the 
electronic  potential  energy  is  a local  functional  of  p(r).9(>  A common  feature  of 
these  models  is  the  replacement  in  (eq.  (5.62) ) of  the  exchange  potential  Vx  by 


Ki  — 4(3/(8jt)p(r)),/J 

where  Vxs  was  introduced  by  Slater93  by  consideration  of  the  properties  of  the 
electron  gas.  A number  of  criteria  have  been  suggested  for  the  choice  of  the  para- 
meter * and  have  been  studied  from  a number  of  viewpoints  including  momentum 


FIGURE  5.5 

Difference  between  the  vx  and  I1F  impulse  Compton  profiles  for  different  choice 
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space  proper ties. -k  ''  l,4  Uf  p was  found  that  1 lie  Compton  profile  was  sensitive  to 
the  choice  of  x In  Fie.  5.5.  the  difference  between  Jx,iq • and  >s  shown1'4  for 
several  choices  of  y including  in  decreasing  magnitude  of  y:  the  original  Slater  value 
a I,  y — 0-8.  the  virial  theorem  value1""  xVt.  and  the  Caspar"  Kohn  Sham1'*’ 
value  of  2 3.  It  is  seen  for  Ar  and  Kr  that  yvl  is  the  best  choice  of  those  illustrated, 
while  for  N'e  there  is  a significant  discrepancy.  Since  the  choice  of  , emphasizes 
\p2>.  it  was  suggested  by  Sabin  and  Smith’*"  that  a different  criterion  he  employed 
for  Compton  profile  studies  w hich  emphasizes  </>  1 >.  namely  ac(.  should  be  chosen 
such  that  JX2lJ0)  = Jtl F(0).  which  is  similar  to  the  method  introduced  by  Bloch  and 
Mendelsohn**9  for  the  choice  of  atomic  shielding  parameters.  In  Table  5.4  we  compare 
impulse  Compton  profiles  and  the  relative  percent  difference  for  y — 2/3,  y.\1  (0-7458), 
and  Ofp  (0-8032).  It  is  seen  that  JXXlr  is  clearly  in  quite  good  agreement  with  ,/HF 
for  small  q and  in  reasonably  good  agreement  for  a larger  range. 

In  addition  to  the  replacement  of  Vx  by  I7*,  there  has  been  much  interest88-91,97 
in  its  replacement  by  various  exchange  and  correlation  potentials  I7xC  based  again 
on  electron  gas  considerations.  We  note  that  several  authors88,95  have  calculated 
l(p)  and  J(q)  from  orbitals  determined  from  different  F*c’s. 


Tattle  5.4  InipiiKe  Compton  profiles"8-95  for  various  local  exchange  potentials  for  neon  ( 1 S).  A gives  the 
difference  from  the  HF  \alue 
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5.4  ELECTRON  C ORRELATION  AND  IMPULSE  COMPTON  PROFILES 

Since  accurate  Hartree  Fock  (HF)  calculations  for  atomic  states  have  become 
available  in  recent  years,  it  is  now  possible  to  test  the  adequacy  of  the  HF  model 
for  the  description  of  electronic  properties.  HF-theory  has  been  found  to  reproduce 
some  expectation  values  of  physical  interest  quite  well,  but  for  others  it  is  not 
sufficiently  accurate.  As  an  example  of  the  latter  category,  one  may  cite  the  electron 
affinity  where  for  atomic  oxygen  the  HF  value  is  of  the  wrong  sign.  The  problem  of 
the  difference  between  the  HF  and  the  exact  non-relativistic  prediction  is  called  the 
correlation  problem.  It  was  originally  formulated  for  the  correlation  energy,'00  namely 
the  difference  between  the  exact  (non-relativistic)  and  the  restricted  Hartree  Fock 
(RHF)  energy,  i.e„ 

£co„  - Fcxiict  - Frhf  (5.64) 

Since  this  error  is  a small  fraction  of  FR,,r.  it  is  usually  only  in  questions  involving 
energy  differences  where  the  correlation  energy  is  important. 
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One  of  the  reasons  for  the  accuracy  of  the  HF  energy  is  expressed  by 
Brillouin's  theorem101  that  F.,n.  and  Ec,  , agree  to  second  order.  Similarly  the 
Moller  l’lesset  theorem  "’’•,0-‘ indicates  that  the  III  expectation  valuesof  one-electron 
operators  should  be  correct  to  second  order  as  well.  1 he  magnitude  of  this  correction 
would  be  dependent  on  the  nature  of  the  operator  and  the  state  involved.  In  addition 
it  is  formulated  for  the  HF  model  without  restrictions  rather  than  the  Rill  model 
which  is  commonly  used. 

Weiss104  was  probably  the  first  to  inquire  about  the  possible  effects  of  electron 
correlation  on  momentum  distributions  and  Compton  profiles.  He  felt  that  atomic 
HI  momentum  densities  would  probably  not  be  as  accurate  as  the  HF  atomic- 
form  factors  which  our  previous  discussion  of  eq.  (5.34)  supports.  Calculations  are 
now  available  at  both  the  correlated  and  RHF  level  for  the  lighter  atoms  (He  Ne). 
In  the  case  of  helium50-45105  the  data  in  Table  5.5  indicates  that  the  111-  and  both 


Table  5.5  Comparison  of  theoretical  impulse  Compton 
profiles  for  helium  ( lS) 


q (a.u.) 

HF30-41 

MCHI'30  105 

Corr15 

00 

1-0705 

1-068 

1 0685 

0-1 

1 0568 

1055 

1 0545 

0-2 

10172 

1-015 

1-0143 

0-3 

0-9557 

0954 

0-9512 

0-4 

08782 

0-876 

0-8737 

0-5 

0-7911 

0-788 

0-7859 

0-6 

0-7004 

0698 

0-6951 

0 7 

0-6112 

0-609 

0-60M3 

08 

0-5271 

0-525 

0-5225 

09 

>4251 

0-449 

0-4465 

1-0 

0-3820 

0-381 

0-3818 

of  the  correlated  profiles  therein  are  in  quite  good  agreement  w ith  a slight  flattening 
of  the  profile  because  of  electron  correlation.  Similar  results  were  found  for  Li 4 
and  other  two-electron  ions.31  Comparison  of  the  HF  and  correlated  / (p)  values 
for  Li  calculated  by  Benesch  and  Smith33  indicates  discrepancies  of  the  order  of 
2 3 per  cent  while  the  Compton  profile  was  more  accurate.  The  maximum  dis- 
crepancies for  J[i/)  occurred  at  the  peak.  2-592S  (RHF)  versus  2-5743  (correlated),  a 
difference  of  0-75  per  cent.  The  atomic  form  factors100  were  in  their  turn  even  more 
accurate,  illustrating  the  comparative  smoothing  nature  of  the  convolution  (eq. 
(5.34) ) and  the  integration  (eq  (5.29))  involved.  The  largest  discrepancies  for  l(p ) 
and  Jiq)  occurred  in  the  valence  or  /.-shell  region  while  in  the  //-shell  region  they 
were  very  small  as  in  the  case  of  Li  * mentioned  above 

Subsequently  Benesch  and  Smith10  and  N'aon  and  CorniIlc,os  examined 
bcivtlium  ( 1 S).  Both  investigations  revealed  very  large  correlation  effects  in  Up) 
and  Jiq)  which  were  primarily  in  the  /.-shell  region.  There  was  a decrease  of 
approximately  7 per  cent  in  the  magnitude  of  J(0|  due  to  electron  correlation.  1 he 
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Table  5.6  Comparison  of  7(0)  values  for  M-n  ral  neon 
I S i wave  fund  ions' 


Wave  function 

JlOl 

£ 

RHF  (dementi) 

2-72?S> 

128  5470 

RMS  (Bagus-Gilberl) 

2-7257 

128  5471 

AH-MCHF 

2 7371 

1 28  S016 

Cl  (VHS)d 

2-7245 

1 2S-K767 

Cl  (YHS-scaled)* 

2-7282 

128  8770 

BCJI' 

2-7340 

to 

* All  calculations  except  BG1  arc  from  Smith  and  Broun  l0<< 

b The  energy  was  not  given  in  Ref  10S  It  is  onlv  stated  that 
VK  per  cent  of  the  correlation  energy  was  recovered  For  com- 
parison. Cl  (VHS|  arid  .AH-MCHF  yield  S9  and  67  per  cent 
respectively. 

‘ The  MC  JIF  w.ive  function  <»f  Ahlrichs  and  Hin/e,,?  uhich 
ir-cluded  onlv  /.  shell  correlation  It  should  he  netted  that  another 
calculation  reported  in  the  literature  * u>ed  an  incorrect  version 
of  tht?  wave  function  1,3 

d The  full  second  order  Cl  of  Viers.  Harris,  and  Schaeffer.1 12 

* Scaled  to  satisfy  the  virial  theorem  in  Ref  109 

' Profile  L.ihul, ifjon  of  Naon  arid  Oornrlle'0*  using  N'esbet*s 
Bethe  Goldstone  increment  method  to  neat  /.  shell  correlation 
only. 


cause  of  the  discrepancy  could  be  identified  as  being  due  to  the  near-degeneraev 
of  the  ls:2s2  and  ls:2/r  configurations. 

The  calculations  on  atomic  neon  by  Smith  and  Brown109  and  Naon  and 
Cornille108  showed  that  for  this  well-closed  shell  system,  electron  correlation  played 
a very  minor  role  in  determining  the  Compton  profile.  This  is  seen  by  comparing 
the  J( 0)  values  which  are  tabulated  in  Table  5.6.  It  should  he  observed  that  all  the 
correlated  values  for  J(0)  with  the  exception  of  one  are  larger  than  the  1 1 F value  and 
hence  provide  a counter-example  to  the  conjecture  that 

*^(0W  > d(0)corr  (5.65) 

This  conjecture  which  is  based  on  the  vi rial  theorem  originated  with  McWeeny  and 
Coulson.38  It  was  suggested  in  recent  years  by  Brown  and  Smith109-110  and  was 
satisfied  by  calculations  for  He.30'4'105  Li.3:  Be,10  los  B.8)  and  H2."p  However, 
recent  calculations  for  H20"  ' which  is  isoelect  ionic  with  neon,  also  contradict  (5.65). 

We  note  that  the  one  calculation  in  Table  5.6  which  obeyed  eq.  (5.65),  did  not 
satisfy  the  virial  theorem  but  when  sealed  to  do  so  no  longer  satisfied  the  conjecture. 
A study  of  atomic  boron  (*/’)  by  Smith.  Brown,  and  Reneseh8 ' revealed  that  correla- 
tion flattened  the  J(<y)  curve  with  the  largest  deviation  (3-23  per  cent)  occurring  for 
J(O).  Here  again,  the  effect  of  electron  cot  relation  on  ./(</)  and  /(/>)  was  largest  in  the 
L- shell  region  and  relatively  minor  for  the  core  (Fig.  5.6). In  order  to  isolate  the 
particular  types  of  electron  correlation  which  are  important  for  describing  / (y?)  and 
•/(<;).  Smith.  Brown,  and  Benesth8’  reported  calculations  for  a scries  of  wave  functions 
which  were  carefully  constructed  by  Seh.wlTer,  Klcmm.  and  Haiiis"4"5  to 
successively  include  the  contributions  in  the  Sinanoglu' 18  decomposition  of  the  exact 
wave  function  into  orthogonal  parts 
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w here  only  one  of  the  substituted  orbitals  belongs  t<>  l!ic  open  shell  I he  remaining 
term  / 1 x , contains  the  remainder  of  the  v.a\e  function  including  the  < Manal 
substitutions  where  all  of  the  substituted  orbitals  are  outside  the  open  shell. 

I'he  Compton  piofiles  coi responding  to  these  correlated  wave  functions  are 
included  in  Table  5 3.  fhe  20  wave  funeiion1 1 consists  of  l1,HIJ1.  - /,M.  i e..  the  two 
configurations  lv  7v7/'  and  Is  ~2/r  w bile  the 'Pol.' wave  function  included  /kHl.  + 
and  the  1-0  (first  order1  ’ 4 or  charge-density  1 "’)  wave  functions  consist  of  <1>KHF  + 
'/.f  h /i\ r • Hie  IJs7  tenn  Cl  includes  most  of  the  same  contributions  as  the  FO 
function  plus  some  of  the  external  contributions.  Inspection  of  the  table  shows  that 
the  latter  two  functions  are  in  \ery  good  agreement  while  the  2C1  and  Pol.  functions 
each  account  for  approximately  half  of  the  discrepancy  at  the  peak.  Smith. 
Brown,  and  Benesch81  concluded  that  inclusion  of  the  internal  orbital  polarization 
and  semi-internal  cot  relation  effects  in  the  valence  shell  is  necessary  to  adequately 
describe  the  Compton  profile  for  open-shell  states  while  the  external  correlations 
play  only  a minor  role.  Since  the  internal  and  semi-internal  contributions  are  absent 
for  well-closed  shell  states,  the  reported  accuracy  of  RHF  for  helium  (!S)  and  neon 
(’S)  is  understood. 

The  most  recent  study  is  that  of  Munch  and  Davidson33  for  vanadium  (4F). 
Inspection  of  Fig.  5.3  show's  that  the  cor  relation  contribution  is  mainly  in  the  valence 
region  w ith  a flattening  of  the  profile  in  agreement  w ith  the  conjecture  of  eq.  (5.65). 


5.5  RF.I.VUMSTIC  PR  Oil  IK  CORK  Ft  HONS 

In  this  section  we  will  discuss  the  effects  of  using  a relativistic  Hart ree  Fock  spinor 
wave  function  on  the  momentum  density  and  impulse  Compton  profile.  We  will  not 
touch  upon  the  important  question  of  a correct  formalism  relating  the  scattering 
cross-section  to  the  impulse  Compton  profile  for  high  incidence  energy  y-ray  photons 
scattered  at  any  angle  except  to  stale  that  it  has  recently  been  reported  by 
Ribberfors1 18  (see  also  chapter  2). 


5.5.1  Relativistic  Hart  ree  Fock  formalism 

Relativ  istic  Hartree  Fock  impulse  profile  calculations  have  been  reported  in  a series 
of  papers  by  Mendelsohn.  Biggs,  and  Mann. ' 911 9,1 20  The  w ave  functions  used  w ere 
Mann's  numerical  relativistic  Hartree  Fock  wave  functions. ,:i  These  were  obtained 
by  solving  the  Dirac  equations  self-consislenlly  taking  all  the  two- particle  inlet  actions 
as  Coulornbic.  In  the  relativistic  central  field  Hartree  Fock  case  a particular  n.l 
orbital  (except  for  s orbitals)  splits  into  two  orbitals,  one  with  j - I + ! and  one  with 
j 1 j.  For  either  of  the  nlj  orbitals,  the  four  component  spinor  may  be  expressed 
in  terms  of  only  two  radial  components:  G„j(r).  a large  component  which  looks  very 
much  like  the  noil-relativistic  radial  wave  function  for  orbitals  of  mu:  II  atoms,  and 
|c|.  a small  component.  In  terms  of  these  components,  the  wave  function 
normalization  condition  may  be  taken  as 
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J i(0'„0( /•))-’  i d.  ! (567) 

The  two  components  of  the  spinet  wh  I - ■ G(r)  l»i| 

harmonics  ol  ordc  I.  1 he  two  components  o!  [he  'j-  ' wl.  : O"  i In  .n. 

multiplied  by  spherical  harmonics  >'f  order  / I lor  the , / l c i>e  1 of « >rd  • 

/ l for  the  j — I - ^ case  Taking  the  Fouriei  transforms  of  the  individual  spinoi 

components,  the  momentum  wave  function  s;. i or  can  d'.>  !■  w.  rti:  . - 1 ; , 

radial  momentum  components  /'•,(;>)  and  /,/_(/>).  where 

■{Zip)  (2  - I 0„ij(r)j\(pr)i  ' dr  (5  ftN) 

J 0 

and 

Xn  jp)  (2  7i 1 1 •’  /•„,;(»■)/.  ir>rh  dr  |5f>‘b 

Jo 

where  i - 1 for  j — l + { and  i — - I for  j l ; and  jApr)  is  a spherical  Bcs-d 
function  of  the  first  kind.  Finally  the  impulse  profile  is  calculated  frot  i 
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1 j \y.nlj\2P'lp 

J n 

i ) MxSA/fl)1  4 {/.ntAp))2',p d/> 


(5  70) 


where  the  factor  F rather  than  2,t.  comes  from  the  normalization  of  G and  F in 
eq.  1 5.67).  A discussion  of  the  numerical  procedures  employed  to  evaluate  t lie  integrals 
may  be  found  in  Ref.  79  (see  section  5.2  1(3)).  In  addition  a.  complete  set  of  re'ativ istic 
Hartree  Fock  Compton  profiles  for  the  individual  orbital-  of  all  free  atoms  with 
Z > 36  is  given  in  this  reference 


5.5.2  Relativistic  corrections  to  the  Compton  profile 

The  most  interesting  result  of  the  relativistic  calculations  is  that  relativ  istic  efTects  on 
outer  electron,  and  therefore  total,  Compton  profiles  are  much  larger  than  one  would 
expect  a priori.  The  orbitals  which  contribute  significantly  at  the  profile  centre  are  the 
outer  electron  atomic  orbitals.  One  does  not  ordinarily  think  of  these  orbitals,  where 
the  electrons  are  moving  comparatively  slowly,  as  being  a fleeted  very  much  by 
relativ  istic  effects.  On  the  other  hand  one  expects  the  relativistic  spatial  wave  functions 
of  rapidly  moving  electrons  in  the  inner  orbitals  of  heavy  atoms  to  be  pulled  in 
considerably  when  compared  to  nun -relativ istic  wave  functions  1 his  pulling  in  ol  the 
relativistic  electron  density  can  be  understood  simply  from  considering  the  m -■ 
velocity  effect  as  a perturbation  to  the  non-rel.ituistic  Hamiltonian  As  such  it  h’- 
the  same  sign  as  the  attractive  nuclear  potential' ' and  the  c'orc  pulls  in  the  c -s.tr  on 
wave  function  so  that  it  is  larger  nearer  the  nude  is.  Since  a more  sharply  pcs  cd 
inner  electron  spatial  wave  function  leads  to  a flatter  momentum  space  wave  lunctiu'  . 


FIG l RF  5.8 

Comparison  of  FLirtree  Fock  [ ) and  relativistic  Hartree  Fock  ( 

momentum  densities  for  the  4s  state  of  lead. 


Clearly  the  relativistic  pulling  in  efTect  is  quite  large.  Another  difference  not  shown 
in  the  figure,  is  that  the  relativistic  electron  density  does  not  exhibit  nodes  since  the 
F and  G components  do  not  go  to  0 at  the  same  point.  A similar  type  of  comparison 
for  the  6 p electron  of  lead  again  shows  a surprisingly  large  relativistic  decrease  in 
J |0)  of  9-52  per  cent,  fn  lead,  as  in  other  materials,  relativistic  flattening  effects  for 
d and  f electrons  are  negligible  for  the  most  part.  When  the  ,/(0|  values  for  all  the 
orbitals  oflead  are  added  up.  one  finds  a relativistic  decrease  of  the  tot  • 1 piofilc  of 
5-24  percent  at  the  profile  centre  Since  these  relativistic  profile  calculation,  were 
performed  within  the  impulse  approximation  using  non  . di/ed  Dirac  spin.  , wave 
functions,  the  normalization  condition  on  the  profile  of  an  . orbital 
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).  ,jt,cd  This  condition  insures  tliat  the  relativistic  . /.I  out  on  tlv  wings  will 

he  1 .1  lh.in  llic  nor-rclativ istic  pr« >lilc  result.  lo  give  the  reader  s >me  estimate 

of  rei.it: i -lie  flattening  results  for  smaller  . loins,  we  cite  the  results  fot  the  rare 
: ises;  for  argon  iZ  IS),  a it:,  tivistic  decrease  in  the  total  J(0|  of  fi  1 S 7 per  cent 
: . f ■ aid:  for  1 ; . pt  t *n  <Z  56|.  a din;  a-e  of  O'.  12  per  cent  is  obtained,  for  xenon 
1/  5 4).  the  eiTect  is  1.47  per  sent. 

5.6  C OMPAUISON  OK  lHhOR\  AM)  I \i  1 KIMKM  FOR  I UKK  Vi  QMS 

X-Ray  Compton  scattering  experiments  using  \1o  Ky.0  7-374  kcV)  and  Ag  Ky  (22- 1 63 
kcV)  sources  have  been  performed  by  F.isenberger  on  helium105  and  neon.'--  y-Ray 
expe:  iiiients  with  a !60  kcV  Fe  source  have  been  made  by  liisenberger  and 
Rs's-d1’5  on  helium,  argon,  and  krypton  Roth  experiments  on  helium  were  found  to 
Vie  in  xerx  good  agreement  with  impulse  Harttee  f ock  profile  calculations,  a not 
vui  pi  i sine  result  since  coi  relation,  binding,  and  relativistic  coi  rections  are  split e small 
for  helium.  The  y-ray  data  for  argon  and  krypton  was  subsequently  re-analysed’  4 
using  a newly  derived  theoretical  relativistic  relationship  between  the  y-ray  cross- 
section  .md  the  impulse  Compton  profile.  Williams1-5  has  shown  that  this  new 
derivation  cannot  be  the  correct  relativistic  one  since  when  it  is  expressed  in  terms 
of  the  usual  v ariables.  it  reduces  essentially  to  the  non  relath  istic  result.  1 bus  w e are 
left  with  a problem  in  interpreting  experimental  profiles  obtained  from  experiments 
using  y-ray  sources,  as  the  experimental  analyses  have  all  used  one  form  or  another 
of  an  incorrect  theoretical  relationship  to  obtain  the  Compton  profile  from  the 
experimental  cross-section.  However,  it  is  expected  that  such  errors  will  not  he 
large.  \s  we  stated  in  the  previous  section,  a new  derivation  of  the  correct  relativistic 
relationship  between  cross-section  and  impulse  Compton  p ofile  has  recently  been 
reported.1"'  Only  at  verv  large  incident  energies  (above  160  ke\  I and  or  small 
scattering  angles  will  the  errors  after  nor realization  aficel  the  profile  to  any  appreci- 
able extent.  In  particular  Ribberfors' 1 h s'jows  that  for  a scattering  angle  of  150 
and  incident  photon  c-c;;.  es  up  t o '50  kc\  . the  corrections  will  no  negligible  out 
to  </  lOa.u  aftei  reno;  mali/ation  It  would  be  most  interesting  to  see  an  analysis 
of  the  experimental  data  using  Ribbci Tors’  iterative  procedure  instead  of  the  tirst 
approximation. 

In  the  rest  of  this  section  theory  and  experiment  for  X-ray  Compton  scattering 
from  neon  will  be  compared  and  some  of  the  apparent  anomalies  which  have  come 
out  of  the  X-ray  experiments  on  krypton  will  be  discussed. 

5.6.]  \ i-i  >n 

In  Table  5.7  s Hisenbergcr's  experimental  JL-sbell  piofilcs  for  MoKa  and  Ag  Ky 
X-rays  scattering  through  1 70  horn  neon  are  eompaied  with  1 H profile  calculations 
and  five  impulse  ca1cuk.‘'ons:  one  UK.  one  suicided  hydiogenic.  and  three  which 
•cl  i.’.e  ia:n’.  boil  i-’T.  .tv  In  the  I'll  c Iculations.  -Tb  ling  factors  Z*(0)  in  the  2s 
and  2/’  wave  f>,  s >.  -.-ic  chi-«cn  to  match  the  1 > , . i - . . 1 se  hydiogenic  (11  1)  piofilcs 
for  the  indiv  dual  states  to  the  icsjvitive  imp. the  llartiee  hock  (lllh)  Complon 
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profiles  at  (/  ■ 0.  With  this  choice  of  effective  charge,  it  was  found  that  the  I ! I /.-shell 
neon  profiles  agreed  with  the  IHF  profiles  to  three  or  more  figures  out  to  a q of 
about  03.  In  argon"'*  a similar  comparison  of  II U with  IH  using  Z*( Of  values 
showed  agreement  of  the  2s  and  Ip  profiles  to  three  significant  figures  or  more  out  to 
a q greater  than  1.  Mendelsohn  and  Bloch’1  asserted  that  if  one  teas  interested  in 
obtaining  an  accurate  FH  profile  calculation  in  the  neighbourhood  of  q 0.  it  was 
most  reasonable  to  use  these  same  7*(  0)  values  which  give  III  profiles  in  ag  cement 
with  IHF  profiles,  not  onlv  at  q 0.  but  also  in  some  neighbourhood  of  the 
profile  centre  This  appears  to  be  the  case  for  FH  Mo  Kv  calculations  for  q between 
0 and  O ',  w hero  the  FT  I results  lie  much  closer  to  the  experimental  profile  than  either 
of  the  two  impulse  calculations.  At  q 0 the  F H J (0)  result  lies  on!}  0-42  per  cent 
above  the  experimental  value.  Since  the  experiment  cites  error  bars  of  I per  cent, 
the  FH  value  falls  well  within  the  error  bars.  On  the  other  hard  the  IHF  7(0) 
value  lies  F93  per  cent  below  the  M oKx  experimental  results  and  therefore  lies 
outside  the  experimental  error  bars.  From  the  FH  calculations  the  value  of  0) 
decreases  as  the  incident  photon  energv  iroc  ises  in  agreement  with  experiment.  1 he 
impulse  profile  calcu'. 'don  has  no  intrinsic  dependence  on  incident  photon  energv 
and  therefore  can  Merer  explain  such  incident  photon  energv  effects  \\  e expect  the 
IHF  to  fall  closer  to  the  \:  Kv  experim ental  resuhs  th  m the  Mo  Kv.  ones  be.  tU'C 
the  IHF  should  improve  in  aeemvcv  with  increased  pit. non  energv.  This  is  indeed 
observed  to  be  the  case  The  \g  •.  EH  /( Ol  value  falls  within  tf  e ; pe  in 
erroi  bars  but  appe  .r  s to  olT-.r  no  real  it.  . t to  IH!  \ pi . > .;h!e  e-  pi  it  ■ m 

of  this  is  given  in  Ref.  19 

We  observe  th  •!  all  four  impulse  A fiat  ms  Ire  whin:  the  ex  pc.  me..  3 er:  nx 
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:il  0 and  that  the  a relatt  -n  cfhvts  on  ./(Of  arc  quite  •>  . -il’  I .r  h ,<{ 

previously  compared  h:s  expci  cntal  results  with  an  ini]  1 atei 

a co! i elated  hast'  <:i  whivh  com  med  ct rors  and  which  y . >e  a ./((*»  .al  ,v 
075  jvr  cent  above  ihe  1111  ic-iilt  I ! is  led  Imn  to  ■ > ’■■  ' ■ ide  that  . ■ \cd 
agreement  with  expsii  .nt  v.ould  he  obtained  when  a better  co: • Fled  -.-..e 
function  vv .is  used.  ! low  c wi  I able  5. 7 - ! < a s that  n -- ; e’a!  ion  effect s on  t lie  ( i e p’nr, 
profile  of  neon  are  t.>o  sit.all  to  be  of  set  voJ  experime  •■tally  at  this  time  or  to  decide 
the  question  of  the  sign  of  the  collection  correction  U>  the  IMJ’  peak  value  As 
Table  5.7  reveals,  the  two  wave  functions  which  yield  the  hugest  correlation  ■ crew 
give  the  opposite  sign  for  the  correlation  correction  to  the  1HF  profile  peak.  As  we 
discussed  in  section  5 3.  virial  the.  rein  scaling  of  the  \ 1 IS  peak  value  raises  it  s-aehtlv 
above  the  IMF  one.  Furthermore,  the  Ml  M(  1 IF  and  BCil  calculations  m Jude 
only  L she!!  correlation  while  the  VMS  function  included  K.I..  and  K I.  ' it. >’,.11 
coi relation.  In  this  ease,  the  experimental  procedure  of  subtracting  the  lllh  core 
contribution  was  followed.  Relativistic  Rattening  of  the  neon  profile  winch  is  not 
included  in  am  of  the  calculations  is  also  about  0.1  per  cent. 

Recently  Bonham  and  associates’  have  performed  experiments  to  determine 
the  Compton  profile  of  neon  fu-m  the  inelastic  scattering  of  75  keV  electrons.  (See 
chapter  8.)  They  find  that  the  1 II  results  lie  closer  than  an)  impulse  calculation 
to  their  experimental  profile,  in  the  vicinity  of  the  profile  centre,  even  though  the  KH 
results  still  appear  to  be  about  3 per  cent  too  low . They  state,1  :i'  ‘the  FH  calculations 
which  explicitly  sum  over  all  transitions  to  the  excited  ion  states  seem  to  present  a 
better  description  of  the  experimental  Compton  process'. 

The  partial  wave  Marlrec  Fock  Slater  (Herman  Skilhnan)  calculations  of 
Grossman  and  Mendelsohn25  on  neon  give  the  closest  agreement  with  experiment  for 
both  small  q (except  for  Mo  Ky  scattering  at  q - 0 where  EH  gives  the  best  result) 
and  large  q for  both  Mo  Ky  and  Ag  Ky.  Compton  scattering.  Of  all  the  theoretical 
calculations,  the  EHS  results  are  found  to  be  in  closest  agreement  with  the  experi- 
mental neon  icsults  over  the  whole  profile. 
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5.6.2  K r»  pf  .1 

In  Tabl  5.8  l impulse  i staff  vistic  Hartree-f  ucli  Mill  t and  the  rel  tivistic 

I H Foek  (1DHF)  calculations  are  compared  with  the  original* 2i  and 

re-an. l!3i  ;-ra\  results  of  I isenhcrg „t  and  Reed  The  only  effects  omitted  from 
the  IDHI  calculation  are  binding  effects  and  correlation  effects.  Cleat l>  bin  1 ug 
effects  will  be  negligible  at  such  high  ir.  idetit  photon  energies  and  large  scattering 
angles.  For  this  closed  shell  system  correlation  effects  are  expected  to  be  quite  small 
Thus  it  appears  that  the  re-  tnu  lysed  experimental  results  are  in  worse  ag  - --.cut 
with  theory  than  the  original  results  near  the  profile  centre 
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